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Abstract 

We consider a special class of solutions of the BKP hierarchy which we call r-functions of hyper- 
geometric type. These are series in Schur Q-functions over partitions, with coefficients parameterised 
by a function of one variable £, where the quantities £(fc), k G Z + , are integrals of motion of the 
BKP hierarchy. We show that this solution is, at the same time, a infinite soliton solution of a 
dual BKP hierarchy, where the variables are now related to BKP higher times. In particular, 
rational solutions of the BKP hierarchy are related to (finite) multi-soliton solution of the dual BKP 
hierarchy. The momenta of the solitons are given by the parts of partitions in the Schur Q-function 
expansion of the r-function of hypergeometric type. We also show that the KdV and the NLS soliton 
r-functions coinside the BKP r-functions of hypergeometric type, evaluated at special point of BKP 
higher time; the variables £ (which are BKP integrals of motions) being related to KdV and NLS 
higher times. 

1 Introduction 

The BKP hierarchy was introduced in [1,2] as a particular reduction of the KP hierarchy of integrable 
equations [1,7]. Like the well-known KP hierarchy, the BKP hierarchy possesses multi-soliton and rational 
solutions. In [3,4], the role of projective Schur functions (Q- functions) [6] in obtaining rational solutions 
of the BKP hierarchy was explained. In fact, the Q-functions are polynomial r-function solutions of the 
BKP hierarchy Hirota equations and these are connected to the rational solutions through a change of 
dependent variables. 

In [9] , certain hypergeometric series in Q-functions (see <|85|) below) were shown to be r-functions of 
the BKP hierarchy. These r-functions are series of the form 

T(t„,tX) = E ^ ZlSx >Q^o)Qx(±t Q *) , (1) 

AeDP 

where £ = {£ m : m = 1, 2, . . . } are arbitrary parameters, £o = 0, Q\ denote projective Schur function, and 
the sum is over the set DP of all partitions A = (Ai, A2, . . . , A&) with distinct parts Ai > A2 > • • • > A^ > 0. 
Considered as a function of the variables t Q = (ti, £3, . . .), series is a BKP r-function, where the set 
t Q are higher BKP times. The second set of parameters t* = (£J,tg, . . . ) give the evolution in a second 
BKP hierarchy. 

Remark 1. Consider a set S, which consists of distinct non- negative integers and includes zero. By DP5 
we denote the subset of all strict partitions whose parts belong to the set S. By the limiting procedure: 
e ?fc _^ iff k does not belong to the set S, we can restrict the sum in |QJ to the set DP5. If S is a finite 
set, then is a polynomial in t Q which describes a rational solution of the BKP hierarchy. 

The typical choice of the BKP higher times is the following: 

JV N' 

mt m = X>3T - (-x k ) m ) , mt* m = J2(yk-(-yk) m )- (2) 
& k 
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In this case the sum in Q ranges over all partitions whose length do not exceed k = min(7V, N* ) . 

We note that a special case of this series QJ , where times were chosen as in |J5}, and e^ m was chosen 
as a step function, was considered in [16] in a combinatorial context, not related to integrable systems. 

In the present paper, we will specialize the variables t* as = (1, 0, 0, . . . ), and study the r-function 
JU as a function of the variables £ m . We find that series Q is a multi-soliton r-function of a different 
integrable hierarchy, which we call the dual BKP hierarchy. The variables £ m of Q turn out to be linear 
combinations of the time variables t Q = (ii , t 3 , . . . ) of the dual BKP hierarchy. We observe that the 
variables £ m (proportional to the times of the dual BKP hierarchy) are integrals of motion of the original 
BKP hierarchy and, simultaneously, the times t Q of the original BKP hierarchy are integrals of motion 
of the dual BKP hierarchy. That is why we call these hierarchies dual to one another. 

The situation we will describe is closely related to corresponding results for the hypergeometric r- 
functions of the KP hierarchy. These r-functions are described as hypergeometric because they generalize 
some known hypergeometric functions of many variables, see [18,19]. We note that the KP hypergeomet- 
ric r-functions yields a perturbative asymptotic expansion for a set of known matrix integrals [10,27]. 
They were also used to construct new solvable matrix integrals [28,29]. Other examples of hypergeomet- 
ric r-functions arise in supersymmetric gauge theories [20], [21], in the problem of counting of Hurwitz 
numbers [22], in counting Gromov-Witten invariants of P 1 [23] and in the computation of intersection 
numbers on Hilbert schemes [24]. In references [30,31], r-functions, which were considered in the con- 
text of c = 1 strings, are also of hypergeometric type. The series for two dimensional QCD, considered 
in [25,26], may be related to the KP hypergeometric r-functions also. We anticipate that applications 
of similar series in Q-functions are found also. 

The series JPl can be studied in the context of random strict partitions. Series generalizes the 
sums over random partitions which are considered in [16]. 

With regard to notation used in this paper, we will use infinite sequences of higher times 

t = (ti, t 2 , t 3 , . . . ), t* = (£*, t^,t%, . . .), (3) 



t = (h,t 3 ,t 5 ,...), t* = (t{,t* 3 ,t* 5 ...), (4) 

and, when they appear as higher times in dual equations will be marked with a tilde. Special cases of 
these, too and t Q (q), will be defined in (|8"0"|) and l|£T|l below. 

2 KP and BKP r-functions 

In this section we will summarise the essential facts about r-functions for the KP and BKP hierarchies 
as given in [1] . The definitions of terms related to symmetric functions may be found in [6] . 

2.1 Schur functions as KP r-functions 

Let A be the complex Clifford algebra generated by the charged free fermions tpi, ip*, where i £ Z with 
anticommutation relations 

fa, ^]+ = [ip*,i>*}+ = 0, fa, ^]+ = S id . (5) 

Consider also the generators 

fcez fcez 

The vacuum expectation value is a linear functional ( ) : A — > C. For linear and quadratic elements in A 
it is defined by (^) = (</>*) = (V^> = (V^*) = and 

m) = l kj i<0 , {^j) = l kj (7) 

x 31 |0 i > J I i < 



2 



For an arbitrary product of linear terms in A, Wick's Theorem gives 



(0\wi ■ ■ ■ w 2n +i\0) = 0, (0|wi • • -w 2n \0) = ^sgn(cr)(0|w CT (i)W CT(2) |0) • • • (0\w a ( 2n -i) w a{2n) |0), (8) 

a 

where Wk are linear terms in A, and a runs over permutations such that cr(l) < cr(2), . . . , a(2n — 1) < 
a(2n), and a(l) < cr(3) < • • • < a(2n - 1). 

The connection between the anticommutation relations and the vacuum expectation value is that 
[wi,ui 2 ]+ = (wiw 2 ) + (w 2 wi). 

For free fermion generators with \p\ ^ \q\, 

WW'(9)) = ^, (9) 
and for higher degree products, Wick's Theorem gives 

W Pl W(qi) ■ ■ ■ HPnW(qn)) = det«Vte)V*fe-)» 

= TT 1 TT ^~ p J^ gi ~ g j) (io) 

r Vi ~ qi i<j ( pi ~ q ^ qi 

Time evolution enters A via the hamiltonian 

oo 

H(t) = J2 H ntn, (11) 



n=l 



where 



For any a e A, we define 



H n = ^kr k+n - (12) 

fcez 

a(t) = e H(t) ae- ff(t) = exp(ad#(t)) a, (13) 



and it may be shown that 

V(p)(t) - expK(p,t))^(p), ^*(«)(t) - exp(-£(g,t))V*(g), (14) 

where e(p,t)=Er=iP^- 

Consider g £ A, which solves the following bilinear equation 



3® 3, 



0, (15) 



where the notation [, ] stands for the commutator, and <g> is the tensor product. Then, one has a r-function 

r(t) = <0(t)>. (16) 

The simplest type of r-functions correspond to multi-soliton solution of the KP hierarchy. Taking 
g = exp(^" =1 aitp(pi)ip*(qi)) gives the n-soliton r-function 



r(t)=det (5ij + -^— e €(p*.*)-«to.t) , (17) 

V Pi 1j ) 

where 

oo 

£(p,t)=5> m t m (18) 

m— 1 

Later, we will also need the soliton solution of the two-dimensional Toda lattice equation (TL) [32], 
which is described by almost the same formula 

r(n, t,t*) = A(t,t*)dct (Sij + _J^_ e £(Pi,".t,n-€fe,™,t,t*) N ) , (19) 

V Pi ~ Qj J 
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where 

oo 

f (p, n, t, t*) = J2 (P™*™ ~ P~ mt *™) + n lQ gP (2°) 

m— 1 

and 

A(t,t*) = e £ ™=i mtm '™ (21) 

which is usually omitted from the definition of the TL r-function [32] since the transformation to nonlinear 
variables removes it from the TL solution. We will also neglect this term for the same reason. It is well- 
known [32], that any TL r-function is a r-function of the pair of KP hierarchies with higher times 
respectively t and t*. There exists a reduction to the one-dimensional TL, which yields also a reduction 
to the nonlinear Schrodinger equation. This reduction is described by the demand that the r-function of 
the two-dimensional TL (up to the irrelevant factor JUJ) depends only on t + 1*. It is provided by the 
condition qi = p7 l in (|19|l and we will use it in what follows. We shall also use the reduction to KdV, 
namely the choice qi — —pi in (|17|l . The KdV r-function depends only on the odd index KP higher 
times, that is, on the sequence t Q . 

Polynomial r- functions are obtained by considering expansions in the parameters pi and qj. First, 
elementary Schur polynomials Sj(t) are defined by 

exp(£(p,t)) = £>(t)p fe . (22) 

fe>0 

Since 

i 

1 =cxp(£(p,t))exp(-£(p,t)) =^^ S! _ 3 (%(-ty, (23) 
we have the orthogonality condition 

i 

J2s*- J (t)s J (-t) = 5^. (24) 

3=0 

For all non-negative integers we can define 

b a 
S(a\b)(t) = (-l) b ^S a +l+fc(t)s fc _ fc (-t) = (-l) h+1 ^ Sfc (t) Sa+b+1 _ fe (-t). (25) 
fc=0 k=0 

This is the Schur function for the partition (a + 1,6-'), which is written using Frobenius notation as 
(a\b). This result is easily proved using the Jacobi-Trudi identity. For any partition function written in 
Frobenius notation, 

S(a 1 a 2 ...a n \b 1 b 2 ...b n )(t) = det (S( ttj \ bj ) ) . (26) 

Using this notation, fHjl gives 

^(t) = ^r Sfc (t)^, #(t) = (27) 

fe>0 k>0 

Consequently, 

i-j 

(v*(t)v;(t)> = J2 ^(t) Se (-t)(^_ k r J+e ) = *k(t)*- 3 --*(-t). (28) 

k,t>0 k=i+l 

Hence we see that 

S(a | b) (t) = (-l) h+1 (^(t)^ h _ 1 (t)), (29) 

that is that S( a |i>) is the KP r-function for g = (—l) b+1 ip a 4>- b -i- More generally, this shows that an 
arbitrary Schur function S(oi— o„|6i— 6 n )i i s a KP r-function for 
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2.2 Q-functions as BKP r-functions 

The subalgebra of A invariant under the symmetry 

^ <-> (-l)ty-i (30) 
is used in a similar way to determine BKP T-functions. There are two bases of neutral free fermions 

0i = ^Wi + (-l)V-i), 0i = - (-l)V-i), (31) 

where igZ, each of which generates this subalgebra. 

Using the results for charged free fermions, the anticommutation relations are 

[<j>i,4>j]+ = [<f>i,<f>j]+ = {-l) % 6i-j, [<j>i,<f>j]+ = 0, (32) 

and, in particular, 0q = (f% = |. Similarly, the vacuum expectation values of quadratic elements are 
given by 

\8j, o t = , (33) 

i > 



(</>;<Aj) = (<j>i<f>j) = < 



and Wick's Theorem is used for arbitrary degree products. 

The neutral free fermion generator is defined by <fi(p) = ^2 neZ P n 4>n- We have (for \p\ ^ \p'\) 

(HpW)) = (34) 



and (<t>(p')<f>(p)) = — (<f)(p)(t>(p')) ■ By Wick's Theorem we get 
(cj>(pi)(j)(p2) ■ --Hpn)) = 



1° 



L V'- ■ ■ (35) 

otherwise 



The connection between the charged and neutral free fermions can be expressed in terms of the 
generators as 

-#(?#*(-?) +ml>{QW{-p) = 0(p)4>(q) + 4>(p)4>(l)- (36) 

In the BKP reduction, even times are set equal to zero and we define t Q = (ti, 0, t 3 , 0, t 5 , . . . ), and 
the hamiltonian 

H B (t )= £ H B t n , (37) 

i>l, odd 

where 

#,f = ^£(-l) i+ W-i-n. (38) 

For the fermion generating function one has 

#p)(t ) - e^C^We-^W = e^^K^^ = e^^fo). (39) 

Note also that 

H(t ) = H B (t )+H B {t ), [H B (t o ),H B (t o )}=0. (40) 
Similar to the KP case, BKP r-functions are defined by 

T B (to) = (h(t )), (41) 

where h is the Clifford algebra of the neutral free fermions (pi. The n-soliton r- function is obtained by 
the choice g = exp(^" = i <H<t>{pi)4>{<li)) ■ 
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The Schur q polynomials are denned by 

exp(2£(p,to)) = $>fc(t )/. (42) 

fc>0 

Thus 

M*o) =^qk(?to)<l>i-k. (43) 

fc>0 

We have 



1 

(&(t )^-(t )) = 2<fc(£t )<&(£t ) + ^(-l) fe % +l (it )^_ fe (it ). (44) 

fc=i 

Since 

1 = exp(2£(p,t o ))exp(-2£(p,t )) = ^©(toJ^-^-t,,) = ^(-^-^(t^-i^y , (45) 
for all n > we have 

n 

^(-i)^i(to)gn-i(t ) = 0. (46) 

4=0 

This is trivial if n is odd and if n = 2m is even then it gives 

m 

q m (t ) 2 + 2 ^(-l) fc g m+fc (t )g m _ fe (t ) = 0. (47) 
k=l 

We can also define 

6 

go,t(to) = q a (t )qb(t ) + 2 ^(-l) fe g +fc( t o)g&-fc(t ). (48) 

k=l 

If follows from the orthogonality condition H46JI that 

q a ,b(to) = -qb,a(t ), (49) 

and in particular, q a>a {t ) = 0. Comparing l|44|) and 1|48[1 . it is clear that 

g«,6(2t ) =2(0„(t o )0 6 (t o )). (50) 

Now consider A = (Ai, A2, ■ • • , A2 n ) where Ai > A2 > •••A 2n -i > A2 n > 0. Note that this is a 
partition with an extra trivial part included if necessary to ensure that the number of parts is even. 
The set of such strict, or distinct part, partitions is denoted DP. For A € DP we define 

Q A (it ) = Pf( gA<>Ai (it )). (51) 

This is the Schur Q-function. By Wick's theorem, 

Q A (±t ) = Pf (2(0 A , (t )<p Xj (t ))) = 2™(0 Al (t o )0 A2 (t ) ■ • ■ cf> X2n (t )>. 

2.3 Hypergeometric r-functions 

These r-functions were introduced by one of the authors in the KP case [8] and the BKP case [9] . 

In the KP case, let r be a function of one variable and for any partition A, let r\(x) = Y[u j) eA r {x — 
i + j), the product being over all vertices in the Young diagram. Then 

r(n,t,t*)=^r A (n)s A (t)s A (t*) (52) 
A 

where t = (ti, t%, £3, . . . ), t* = (t*, £3, £3, . . . ), is the KP hypergeometric r- function. 



G 



In the BKP case, r\ has a different definition: if A = (Ai, . . . , Afc) then r\ — ]X=i r (l) r (2) ■ ■ 1 
If we introduce new variables r(k) = e^ k ~^ k - 1 , £_i = then r\ = Yll=i e ^ Ai 
With these definitions, 

T(t ,tx)= E 2-^v A Q A (|t )Q A (it:) = J2 2- £(A) n eCA '^(5 t o)QA(K) (53) 

AeDP AeDP i=l 

is the BKP hypergeometric r-function. By Remark^one can restrict sum (|53|l to the sum which ranges 
over DP 5. 

It can be show that this is a r-function since 

T(t ,£X) = (0|e E ^i.odd^t» e E~ = _«, (-) n €»:<A„<A-„: e E i > 1 ,odd»-nC| ) (54) 

3 Main results 

3.1 KP infinite soliton solution 

Let 

3 sol = exp I OiJ^WiQi) 1 ( 55 ) 



\i,j>0 



Then 



^_sol . /„sol 



5 sol (t)) (56) 
1 + ^2a i::j (ip(p l )ip*(q j )(t)) + ^2 ( a i,j a k,l - <H,iaj,k){ip(Pi)il'* (q 3 )^(PkW («)(*)) + ' ' ' 



A 1 



. e «(pi,t)-5( 9j ,t) 



3 J Pi - 1i 



\^ i (Pi-Pk){qj -gQ c g(p i ,t)-g(q 1 -,t)+g( Pfc ,0-g(q i ,t) + ... 

^ * fe) v (Pi ~ ?j)(Pfc - ft) (P« - - Pfc) 



i,j,k,l 



where 



(57) 



denotes the 2x2 minor of the infinite matrix (ojj) containing the ith and jth rows and the kth and Zth 
columns. 

If a,i t j = S(j|j)(t*) then the coefficients can be written as 

sx(t*) (58) 

for all partitions A. 

3.2 KdV soliton solution 

The KdV reduction of the KP soliton solution l|17|l , 

ffKdv = ex P ^2a>iPii>(PiW(-Pi) , (59) 

\i>0 I 
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gives rise to the following soliton r-function, which depends only on the higher KP times with odd 
numbers t Q = (ii, t 3 , t 5 , . . . ), 



rk% := (flSSv(t.)) = det + ^ e ^M~^M 

V Pi + Pj 



(60) 



1 (Pi-pj) 2 cwt+m 

(p. + P.J') 2 

1 (g» - Pj) 2 (p» " Pfc) 2 bj - Pfc) 2 c?n +„,+„ fc 



where 



Vi 



2J2 t 2m -iPi m ~ 1 + log ai , 1 = 0,1,2,... (61) 



m— 1 



Remark 2. We note that the fractional linear transformation of the (complex) plane of spectral param- 
eters pi e C 

api + b ad -be . 

— , —— = ±1, 1 = 1,2,... (62) 

cpi + a ad + be 

leaves invariant the factors 

{ PZ P € , i,i = l,2,.... (63) 
It is easy to see that there are two distinct types of invariance; a = d = 0, b, c =/= and a, d ^= 0, b = c = 0. 

3.3 NLS and one-dimensional Toda lattice soliton r-function 

Let us consider the reduction of the TL soliton r-function (|19|) to the one-dimensional Toda lattice 
(1DTL) reduction, which is qt = p" 1 , see [32]. If, in addition, = 1, then, it is also a reduction to the 
nonlinear Schrodinger equation (NLS). For the multi-soliton tau function we have 



Sidtl = exp yfe -ft 1 )^ W(pr X )j 



(64) 



7fg m (n,t,t*) = 7$ix(n,t + t*) := 

= (rfgn.M.t*)) = det (s itj + ^-^ ^MW^n] 

_ i + V ie 1 " 4- V 1 &-Pj) 2 e m+vi 
l + 2^2 e + 2^22 {piP3 - 1) 2 e 

, \- 1 (p l -Pj) 2 (p t -Pk) 2 (Pj~Pk) 2 m +r lj +r lk ffifl 
+ ^ fc 23(^-l)2(^-l)2(^_l)2 e +' > ^ 

where t = (ii, i 2 ,t 3 , . . . ), t* = (if, f|, t|, . . . ), and 

OO 

'7,= ^(pr-pr^^ + O + anlogft + bgai , i = 0,1,2,... (66) 

m— 1 

For the nonlinear Schrodinger equation the n-dependence of the r-function is irrelevant. 
Remark 3. Here, the fractional linear transformations 

j_ aPi + b . 

Pt~*±7 ; — , « = 0,1,2,... (67 

bpi + a 

where a and b are not both zero, leave invariant the factors 

fc^C , i,i = 0,l,2 (68) 

(PiPj - l) 2 
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3.4 BKP infinite soliton solution 

Now writing qi = —pi and choosing skew-symmetric matrix entries dj% — —a>ij in (|55(l gives 

5 so1 = exp Hj(-P^(PiWi-Pj)+P^(PjWi-Pi)) ■ (69) 

\0<i<j J 

By H36|l this may be rewritten as 

g so1 = exp ( Y a ^ (^(Pi)^(Pi) + 4>(Pi)kPi)) ■ (70) 
\0<i<3 ) 

Since [(pi, 4>j]+ = 0, [4>i4>k, <j>j<j>i] = and so we can factorize as g = hh where 

h = exp I Y Hj^tPiWiPj) J > ^ = CX P I X] a i,3^(Pi)^(Pj) ■ ( 71 ) 

\*<J / \*<J / 

Then, we have 

r| ol (t ) : = (h(t )) = (h(t )) 

= 1+ Y °*^^(P*)^(Pj)(to)) 
0<i<j 

+ X! ( a i,i a k,i ~ ai,kOj,l + a%,iaj,k){<l>(Pi)<f>(Pj)<l>(pk)<f>(pi){to)) H 

0<i<j<fe</ 

l Pi --Pj r £( Pi ,t„1+£fa-,-,tJ 
2pi+ pj 



0<i<j 



o<i^*<i 22fe+ft)(p fc +p0 



where 



(i,i,*,0» (73) 
denotes the pfafhan minor of the infinite skew-symmetric matrix (a^ ) containing the ith, jth, fcth and 
Zth lines. Using gives 

r sol (t )=r| ol (t ) 2 . (74) 
If aij — <7i.j(it*) then the coefficients can be written as 



Qxm) (75) 



for all partitions A into distinct parts. 
Remark 4. The factors 



Pi+pj & + Pj){pk + pi) 

in H72|) . are invariant under the transformation 

api + b 



cp t + d 



i = 0,1,2,..., (76) 



atj -> 7aj,j , i, j = 0, 1, 2, . . . , (77) 
where a, 6, c, d are any complex numbers such that 

ad + be 



ad — be 

is not zero or infinity. 



(78) 
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Remark 5. Although the free fermion generator is not defined if its parameter is 0, i.e. </>(0) does not 
make sense, the limit 

lim MW)) = \ (79) 

as given by 135(1 . is well defined. 



3.5 Useful Lemma 

Let us introduce the following notation: 

too = (1,0,0,0,...) , (80) 

and 

2 

t {q) = (h(q),t3(q),t 5 (q),...) , t 2m -i(q) = _ jwj- _ g2m _^ , <2m = 0, m=l,2, ... (81) 

Remark 6. Let us notice that t^ can be viewed as given by (J2J), where we take x\ — x-i = ■ ■ ■ = xn = N^ 1 
and N — > oo. Similarly, t Q (g) is given by 10), where Xk — q k ~ 1 ,k = 1,2, ... . As for too, if / satisfies 
/(cii,c 3 t 3 ,c 5 i s , . . . ) = c d f(t 1 ,t 3 ,t 5 , . . . ) for some d G Z, we have h d f(t (q)) -> /(too) as ft := logg -> 0. 
In that sense 1 (80(1 may be considered as a limit of 1(81(1 as q — > 1. 

We have 

Lemma 1. Let A = (Ai, . . . , A&) be a strict partition. Then 

«»*.>- n £ n£^. (82) 

i— 1 i<j J 

and 

q x * - q Xj i r (-g; g) 



where 



(p; g) m := (1 - p)(l - pg) • • • (1 - pq^ 1 ) (84) 



3.6 Hypergeometric functions related to the projective Schur functions 

We consider hypergeometric r- functions JD, where we specialize the variable t* respectively by lf%U|l and 

r{toAM= E 1^11^' (85) 

AGDP S A i=l 

r(t 0> e,t (g))= £ |^n^S (86) 

where 

#a = 0A(it oo )~ 1 , flJ( g ) = Q A (§t (g)) _1 (87) 

are respectively so-called product-of-shifted-hook- length [6], which generalize the notion of the factorial 
for strict partitions and its g-analog (shifted hook polynomial) . We took into account Remark ^ to 
restrict sums over all strict partitions to the subset DP5. 

In the case that DPs is the set of all strict partitions, namely, DP, the series 185(1 and 1(86(1 may be 
considered as multi- variable generalization of hypergeometric function (respectively basic hypergeometric 
function), which we obtain when £(X) — 1 and t D is of form (J2J where Ni = 1. 
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The notation Q(x^ N ') below will be used for Q(^), where t Q is defined by J2J. Let all parameters bu 
be not equal to negative integers. Let in (|85|l we choose 



Then tau function (|85|l defines the following hypergeometric function 



(88) 



xl^P llfc=l(°fe)A W A 

i (A<N 

which generalizes the hypergeometric function of one variable 

p F s (a u ...,a p ;h,..., b s ; * (1) ) = £ ^ ?* " ^ (90) 

The function Ij89(l was introduced in [9]. Here we introduce the g-deformed version of (|89J) . If in (|8fj|> 

we choose 

f n _ nti(g a -;g)n (Q ) 

ni=i(^;?) n ' 1 j 

tau function l|86(l defines the hypergeometric function 

$ r ffl1 o-6i V 2 ^(A) nLi(g afc ;g)AQA(xW) 

which generalizes the basic hypergeometric function of one variable 

°° TIP („0. 



Fin a-b b-rW) V M^llk ^ 



(93) 



Hypergeometric functions (|%9")l and may be also considered as multisoliton tau functions, see next 
subsection. 

3.7 Soliton solutions and rational solutions 

Let us recall, that in Remark^Jthe set of the partitions DPs was defined via a set of distinct non-negative 
integers S which includes zero. 

Theorem 1. Let r(t , £, too) be defined by l|85jl ■ and r sol (t, t*) 6e defined by l)72[l. where 

am + b 

Pm = —5 , ™eS (94) 

cm + a 

(m particular, one can take integer momentum p m = m), with any a, 6, c, d swc/i that 

ad + be 
ad — be 

is not zero or infinity, and 



1 , (95) 



^•=7^(3*0) (96) 
Let for a given set of the numbers p m , m € S, the variables t,t* are related to the variables £ as 



k=l 

Then we have 



^2(PrJk-p m %)+logml (97) 
=1 

r sol (t) =r(t ,£,too). (98) 
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Proof. Let us compare (|85|l and (|72|) . First replace t Q with t D in 1|72|) . Then a typical term on the right 
hand side is 

where < rt\ < n<i < ■ ■ ■ < and k is even. We can set £(pfc,t ) = £fc, choose the parameters p^ = k 
and the pfaffian elements to be (for i < j) 

2gj(|t Q ) . = Q 

(/./) = -( rt /i, s , (ioo) 



so that 



nti ^ QA(|t ) »i > o 



(ni,n 2 , . . . ,n fc ) - fe i , ljL N ^„ (101) 



where A is the partition into distinct parts (nk,nk-i, ■ ■ ■ ,ni). 
Thus the typical term may be written as 

2 ^(A)-Q e « fc Q A (i to )g A (i tco)j (102 ) 

i=l 

for any partition into distinct parts. The partitions into an odd number of distinct parts come from those 
terms for which n\ =0. The extra factors of 2 in the pfaffian element in \ 100(1 and l(101[) are needed 
because in the case ri\ — 0, the length of the partition £(\) is A; — 1 not k. 

This establishes the connection between l|85[) and 172fl . □ 

The hypergeometric function (|89|l is an example of multi-soliton tau function of the dual BKP hier- 
archy, evaluated at special values of times t„, see (|88|l and (|97(l . 



Theorem 2. Tau functions T(t ,^,t 00 ) and Tjgj V (t ) are defined respectively by (|85|l and by (|60() . Let 



-sol 

una i yi>o: C,; L oo ) wmL r 

us choose p m in ()60JI by 



am + b r- Q / 1 nQ^ 

Pm = — -j , me S , (103) 

cm + a 



where a, b, c, d are any parameters with the property 

ad — be 



ad + be 



= ±1 (104) 



(in particular, one can choose integer momentum p m — m, m £ S). Let the numbers £ m in l#,5)) be 
related to i](t 0l p m ) in (|60|) by 

oo 

£ m - log m! = T) m := 2 J2 Pm" 1 ^-! (105) 

k=l 

Then 

^Kdv(to) = r(t 00 ,e,t 00 ). (106) 

The hypergeometric function l|89|) . where x\ = Xi = ■■■ = N^ 1 and N — ► oo, is an example of 
multi-soliton KdV tau function, evaluated at special values of times t Q , , see (18811 . (|105ll and Remark H3 

Theorem 3. Tau functions r(t , £, t (g)) and T^ TL (h, t, t*) are defined respectively by (|86(l and 
(|65|l . Let p m in 165fl oe chosen by 

p m = ±^±R, m eS (107) 
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(in particular, p m = q m , m G S). Let the numbers £ m in |#6|) be related to ?/(t ,p m ) in i|65|) by 

U - log / g;g) 7 = ffc, := 2 f>* - p-*)(t fc + tj) + 2n logp m (108) 

I <7,<7V fe=1 

Then 

rfSWM + t*) = r(t (g),£,t (g)). (109) 
-For i/ie particular choice p m = m, m£ 5, i/ws is ako -/ViS 1 multi-soliton tau function. 

The hypergeometric function l|92|l . where = q k ~ 1 , k = 1, 2, . . . and -/V — > oo, is an example of multi- 
soliton 1DTL tau function evaluated at special values of times t + 1*, see H91(l . I|108(l and Remark 

Remark 7. In case S is a finite set, the polynomial r-function of type (|85|l (and iJHSJ) is related to the 
soliton r-function with a finite number of solitons. 

Remark 8. We note that the higher times t D of the BKP hierarchy we started with are integrals of motion 
for (solitonic) r-function (|72|l of the second BTL hierarchy. Simultaneously, the higher times n, t,t* play 
the role of integrals of motion for the original BKP hierarchy. We therefore call these hierarchies dual 
to each other. 

Remark 9. An oo-soliton solution with spectral parameters lying on a lattice appeared in [12-14] in a 
different way and in a different context. Other links between soliton and rational solutions of the KP 
hierarchy were found in [15]. 



4 Discussion 

An interesting problem is to study the asymptotic behaviour of hypergeometric r-functions. We hope 
to apply methods of soliton theory to conduct this study. We note that the asymptotic behaviour of 
infinite soliton r-functions, similar to those considered in the present paper, was studied in [13]. 

We hope to apply the series to certain problems. 

(1) Let us consider an integral 

'<».<-*> -wU-nn t4U'4 R (110) 

where T is a integration domain in each (zk, z£) plane (k = 1, . . . , N), and /i is a function such that 

fi(zz*)z n dzdz* = I I [i(zz*)z* n dzdz* = 5 n fi, (111) 

and 

r r 

fi(zz*)z n z* m dzdz* = 2S n<m e in . (112) 



The series (in the case that the sum ranges over partitions of length £(X) < N) yields the 
asymptotic expansion for the integral H110(l [9] 



7(JV,t OJ t;)= J2 2-'W e E«ii^Q A (|t )Q A (|t:). (113) 

A£DP,Z(A)<JV 

The restriction £(X) < N makes the difference between the r.h.s. of (I113|l and r(t Q ,£,t*). 

In the limit N ~ * oo (which is typical for applications), the restriction £(X) < N is irrelevant for the 
perturbation series related to (f 1 1 ()|> . and therefore, this series coincides with (JJJ. Also, in the case that 
at least one of the sets t Q , t* has the form of 

JV JV 
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then by the bosonization formulae and Wick's Theorem we get that 

Q A (it o ) = 0, £(X) > N (115) 

Therefore, in this case, the integral I(N, t OJ t*) is the BKP r-function r(t ,£,t*). 

It may be interesting to apply to matrix models and to statistical models where partition functions 
reduce to integrals (|110|l . see, for instance, [11] for examples of integrals similar to (|110(l . 

It is interesting to compare integrals (|110|) with supersymmetric matrix models [33] . 

(2) The series JQl can be studied also in the context of random (strict) partitions. 

Random strict partitions were considered in [16] and, in particular, the "shifted" measure Q \(x)Q \(y) 
on (strict) partitions, were considered in [16]. In this paper, the series where all £ rl = 0, and where 
Ai does not exceed a certain given number was studied. 

Let us remark that the expression 

n«* (n^n^) = {-k)'^- (116 > 

(where is known to be related to the number of shifted tableau of the shape A, see [17], [6]) in the 
case 

£„ = 0, n = 0,1,2,... (117) 
may be considered as the analogue of the Plancherel measure [34] , while in the case 

e f » = (z) n (l - Z) n , (z) n = r( ^, + " } , 71 = 0,1,2,... (118) 

T(z) 

as an analog of the so-called (z)-measure on random partitions (see [34]). 

(3) Finally, we note that KdV soliton solutions with integer momenta was first considered in [12] . Each 
KdV solution of this type yields a wave operator df — A 2rl _i — u, where the potential u — 2dii T Kdv(h)> 
which satisfies the so-called generalized Huygens principle. 

Acknowledgements 

The authors are grateful to T. Shiota for useful discussion. The work was supported in part by the 
Russian Foundation for Fundamental Researches (Grant No 02-02-17382), and the Program of Russian 
Academy of Science "Mathematical Methods in Nonlinear Dynamics" . We thank Claire Gilson and Chris 
Athorne for the organization of the visit of one of the authors (A.O.) to the University of Glasgow, during 
which work on this paper was initiated. 

References 

[1] E. Date, M. Jimbo, M. Kashiwara and T. Miwa, "Transformation groups for soliton equations", in 
Nonlinear integrable systems — classical theory and quantum theory eds. M. Jimbo, and T. Miwa, 
World Scientific, pp. 39-120, (1983). 

[2] E. Date, M. Jimbo, M. Kashiwara and T. Miwa, Physica 4D, 343, (1982). 

[3] Y. You, "Polynomial solutions of the BKP hierarchy and projective representations of symmetric 
groups", in Infinite- dimensional Lie algebras and groups, pp. 449-464, Adv. Ser. Math. Phys., 7. 
World Science Publishing, Teaneck, New Jersey, (1990). 

[4] J.J.C. Nimmo, J. Physics A, 23, 751, (1990). 

[5] I. Schur, J.Reine Angew. Math. 139, 155, (1911). 

[6] I.G. Macdonald, Symmetric Functions and Hall Polynomials, Clarendon Press, Oxford, (1995). 

[7] V.E. Zakharov and A.B. Shabat, J. Funct. Anal. Appl. 8 226, (1974). 



14 



[8] A.Yu. Orlov, "Hypergeometric r-functions r(t,T, t*) as oo-soliton T-function in T variables", non- 
tin. SI/ 0305001. 

A.Yu. Orlov, Theoretical and Mathematical Physics, 137, 1573-1588, (2003). 

A.Yu. Orlov, "Soliton theory, Symmetric Functions and Matrix Integrals" , SI/0207030 vl (To appear 
in Acta Applicandae Mathematicae, 2004). 



I.K. Kostov: "Solvable statistical models on a random lattice" hep-th/9509124 v4. 

Yu.Yu. Berest and I.M. Loutsenko, Comm. Math. Phys., 190, 113, (1997). 

I.M. Loutsenko and V. Spiridonov, Journal of statistical physics, 99, 751, (2000). 

I. Loutsenko and V. Spiridonov, "Self-similarity in Spectral Problems and q-special Functions" , 
Nuclear Phys B 538 (1998) 731, solv-int/9909022 



A. Mironov, "Seiberg-Witten theory and duality in integrable systems", hep-th/0011093 



C. Tracy and H. Widom, "A Limit Theorem for Shifted Schur Measure" math.PR/0210255. 

J.R. Stembridge, J. Alg. Combinatorics, 1, 71-95, (1992). 

K.I. Gross and D.S. Richards: Trans. Amer. Math. Soc, 301, 781, (1987). 

S.C. Milne, "Summation theorems for basic hypergeometric series of Schur function argument", in 
Gonchar, A. A. and Saff, E.B. (eds) Progress in Approximation Theory, pp. 51-77. Springer- Verlag, 
New- York, (1992). 

N.A. Nekrasov: "Seiberg-Witten Prepotential from Instanton Counting", hep-th/0206161 vl. 

A.S. Losev, A. Marshakov and N.A. Nekrasov, "Small Instantons, Little Strings and Free Fermions", 
hep-th/ '0302191 v3. 

A. Okounkov, "Toda equations for Hurwitz numbers", AG/00044128 vl. 

A. Okounkov and R. Pandharipande, "The equivariant Gromov-Witten theory on P 1 ", AG/0207233 
vl. 

Wei-Ping Li, Zhenbo Qin and Weiqiang Wang, "Hilbert Schemes, Integrable Hierarchies and 
Gromov-Witten Theory", AG/0302211 vl. 

A. Migdal, Zh. Eksp. Teor. Fiz., 69, 810, (1975). 

K. Saaidi and M. Khorrami, "Large N limit of the nonlocal 2D Yang-Mills and generalized Yang- 
Mills theories on a cylinder" , [hep-th/ 0304204 \ vl. 

J. Harnad and A.Yu. Orlov, Theoretical and Mathematical Physics 137, 1676-1690, (2003). 

A.Yu. Orlov, "New Solvable Matrix Integrals", nlin. SI/0209063 (To appear in Int. J. Math. Phys. A 
in 2004). 

A.Yu. Orlov, "New Solvable Matrix Models — U(n) case" in Nonlinear Physics: Theory and Exper- 
iment. II, pp. 99-100, Universita' di Lecce - Consortium Einstein; editors: M.J.Ablowitz, M.Boiti, 
F.Pempinelli, B.Prinari; World Scientific New Jersey, London, Singapore, Hong Kong (2002). 

T. Nakatsu, K. Takasaki and S. Tsujimaru, Nucl. Phys. B, 443, 155, (1995). 

K. Takasaki, Comm. Math. Phys. 181, 131, (1996). 

K. Ueno and K. Takasaki, Adv. Stud. Pure Math. 4, 1-95, (1984). 

T. Guhr, J. Math. Phys. 32 (2), Febraury 1991 

A. Borodin and G. Olshanski, "Random partitions and Gamma kernel", math-ph/0305043 



15 



